Temperature probes in binary granular gases 
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We investigate the validity of Fluctuation-Dissipation (FD) relations for a mixture of two granular 
gases with different physical properties (restitution coefficients or masses) subject to stochastic 
driving. It is well known that the partial granular temperatures T\ and T2 of the two components 
are different, i.e. energy equipartition is broken. We observe, with numerical simulations of inelastic 
hard disks in homogeneous and non-homogeneous situations, that the classical equilibrium Green- 
Kubo relations are satisfied separately for each component of the gas, the role of the equilibrium 
temperature being played by the granular temperature of each component. Of particular interest 
is the limit in which one of the two components consists of only one particle, representing a non- 
perturbing thermometer. In this case it turns out that such a thermometer is measuring its own 
temperature and not that of the surrounding granular media, which in general will be different. 
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INTRODUCTION 

Granular gases pj, i.e. gases of particles interact- 
ing through dissipative collisions, represent an important 
paradigm for the study of non-equilibrium stationary 
states. Due to the dissipative nature of the interactions, 
granular gases have to be considered as open systems 
and therefore concepts from equilibrium thermodynam- 
ics cannot be applied, at least in a straightforward way. 
However, by analogy with molecular gases, a "granular 
temperature" T g can be denned in terms of the kinetic 
energy per particle. If the system is driven by an external 
energy input a stationary state is reached: T g fluctuates 
around a given value, while a flux of energy from the driv- 
ing source goes into the system and is dissipated through 
inelastic collisions. Experiments, theory and simulations 
have recently shown that in a mixture of different grains, 
as soon as the inter-particle collisions dissipate energy, 
the equipartition of energy among the mixture compo- 
nents is lost Ell S H MM II 11 El 13 113 even in the 
tracer limit [13j. Although in sharp contrast with the 
behaviour of molecular gases at equilibrium, this viola- 
tion is not surprising because in a generic open system 
the equipartition of energy is not expected. However, it 
is natural to ask whether T g could have some different 
"equilibrium" meaning, or if it is just a measure of veloc- 
ity fluctuations. A hallmark of equilibrium phenomena 
is the well-known fluctuation dissipation theorem, relat- 
ing the response of a system to a perturbation to the 
corresponding correlation function measured in the un- 
perturbed system: the response of an observable B at 
time t to an impulsive perturbation Ha at time t = can 
be obtained as 



5{B) 

SllA 



ld_ 



■B(t)A(0)) 



(1) 



where A is the observable conjugated to \ia and T is the 
equilibrium temperature of the system (brackets denote 
averaging over thermal history). A recent investigation 
has shown that, for non-equilibrium monodisperse driven 
granular gases, this relation is still obeyed if the equilib- 
rium temperature is replaced by the granular tempera- 
ture 0. 

A natural question then arises: what happens to the 
relation Q in a mixture which displays more than one 
granular temperature ? In this Letter, we will address 
this question in the case of a binary mixture subject to 
a homogeneous driving. It turns out that each species 
obeys an FDT relation with its own granular tempera- 
ture |Fo| as proportionality factor. This result is partic- 
ularly relevant to experiments where a tracer, immersed 
in the granular gas and acting as a probe (and in par- 
ticular as a thermometer), may have different phy sical 
properties with respect to the surrounding gas |17| . In 
equilibrium measurements, the temperature does not de- 
pend on the thermometer. For granular materials, in- 
stead, the thermometer measures its own temperature, 
in general different from that of the surrounding gas. 



MODEL 

We consider a volume V containing a mixture of N = 
N% + N2 inelastic hard spheres (IHS) in dimension d = 2, 
Nx and N 2 being the number of particles in component 1 
and 2 of the mixture, respectively. The spheres have di- 
ameters a (identical for the two species) and masses m Si 
(where 1 < i < N and Si is the species index, I or 2, of 
particle i). In a collision between spheres i and j, char- 
acterized by the inelasticity parameter called coefficient 
of normal restitution a SiSj , the pre-collisional velocity of 
particle i, Vi, is transformed into the post-collisional ve- 
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locity v[ such that 

v[ = Vi ^ (1 + a StS )(a ■ v tj )a (2) 

m st + m Sj 

where Vij = — Vj and <r is the center to center unit 
vector from particle i to j (a SiSj = a SjSi so that the 
total linear momentum m Si Vi + m Sj Vj is conserved). The 
granular temperature of species s is given by its mean 
kinetic energy T s — (m s i)g)/d. 

The loss of energy due to collisions can be compen- 
sated in various ways. In experiments the energy is 
typically supplied at the boundaries, leading the sys - 
tcm to a heterogeneous stationary state 0, 0, IS Ell • 
In order to avoid the complication of strong tempera- 
ture heterogeneities, we will use a homogeneous driv- 
ing in the form of a "thermostat" : in this mechanism 
(which recently has attracted the attention of many the- 
orists [H US IH IH ES HI), the particles are sub- 
mitted, between collisions, to a random force in the form 
of an uncorrelated white noise (e.g. Gaussian) with the 
possible addition of a viscous term. The equation of mo- 
tion for a particle is then 

mi— = Fi + rriiRi - QVi (3) 

where Fi is the force due to inelastic collisions, £ is the 
viscosity coefficient and (Ri a (t)Rjp(t')) = ^o8ij8 a p8{t — 
where Greek indexes refer to Cartesian coordinates 

At the level of Boltzmann kinetic equation, the tem- 
perature ratio of a binary granular mixture subject to 
stochastic driving of the form given above has been ob- 
tained in 

7] 

for the case Q — and in [HI for C + 0. 
METHODS 

We have used two different simulation methods: the 
Direct Simulation Monte Carlo (DSMC) |H |H which 
neglects pre-collisional correlations and therefore enforces 
the molecular chaos hypothesis (factorization of the two- 
particles distribution functions) and Molecular Dynam- 
ics (MD) simulations. The DSMC algorithm can be used 
to simulate the homogeneous Enskog-Boltzmann equa- 
tion, i.e. any particle can collide with any other in 
the whole volume V, or to simulate an a priori non- 
homogeneous system, by dividing the space into small 
cells (of linear size smaller than the mean free path) and 
enforcing molecular chaos in every single cell. On the 
other hand, molecular dynamics (MD) simulations inte- 
grate the equations of motion of the IHS: in this case 
we consider N hard spheres in a square box of linear size 
L = VV, with periodic boundary conditions, random ini- 
tial velocities, and we use an event-driven algorithm to 
study their dynamics. All the investigation methods used 



have shown similar results. We have used sizes N = 500 
to N = 5000, and averaged over up to 10000 realizations 
of the dynamics. 

Our analysis relies on two sets of independent measure- 
ments, i.e. two choices of the pair response-correlation 
to be measured. The first one consists in a measure of 
mobility and diffusivity. The mean-square displacement 
(MSD) for each species, 

N l(2) 

B m (t,t') = —- £<|r J -(t)-r j (f)| a ) , (4) 

behaves as ~ 4Z? 1 ( 2 )(£ — t') for large time differences. 
The mobility of a tracer particle is measured by apply- 
ing a constant and small (2^ drag force £ = £e x to a 
given particle, labeled 0, for times t > t'. Due to the vis- 
cosity induced by the collisions with other particles, the 
perturbed particle will reach at large times a constant 
velocity /x, related to the response by 

Xi,2(t,t') = -((r (t)-r (f))-e x )^ f it at large t. (5) 

By successively using as test particle one particle of each 
species, one obtains the two responses xi and %2, and 
thus the mobilities fi± and /i2- Two Einstein relations 
(/is = 2Di/Ti) can then be checked, e.g. by plotting Xi 
vs. Bi. 

Another totally independent way of checking FD rela- 
tions in granular gases has also been used in 14]: once 
a steady-state has been reached, the system is perturbed 
impulsively at a given time to by a non-conservative force 
applied (non-uniformly) on every particle (we will take 
<o = without loss of generality). The response is then 
monitored at later times. The force acting on particle i 
is 

F(r,,t)= 7i €(r i ,t) (6) 

with the properties Vx(/0,V^ = 0, where 7j is a 
particle dependent variable with randomly assigned ±1 
values. A simple case is realized by a transverse per- 
turbation £(r, t) — (0, £ cos(k x x)S(t — t )) 27], where k x 
is compatible with the periodic boundary conditions, i.e. 
k x = 2ixnkjL x with integer and L x the linear horizon- 
tal box size. The staggered response function (i.e. the 
current induced at t by the perturbation at to), and the 
conjugated correlation, 

R(t,t ) = liVi(t) cos(k x Xi(t))) , 

' i 

C(Mo) = (^2yi{t)y i (.t )cos{k x {x l (t) - x l {t )}}) 

i 

are related, at equilibrium, to the correlation by the FD 
relation R(t,to) = ^C(t,t ), T = 1/(3 being the equilib- 
rium temperature. 
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Ref. 01 has shown the validity of this relation in the 
context of a monodisperse granular media heated by a 
thermal bath with temperature Tf,, reaching in this way 
a stationary state with granular temperature T g < T^. 
In this case the FD relation holds by replacing T with 
T g . For a binary mixture, two sets (Ci,Ri) (i = 1,2) of 
correlation and response are measured separately, thus 
obtaining two plots Cj vs. Ri. 

There are different reasons to test various pairs of re- 
sponse and correlation functions. First, at equilibrium 
FDT links any couple of conjugated correlation and re- 
sponse function with the same temperature. Out of equi- 
librium, it is possible, a priori, that FDT could be valid 
for some observables and not for others. Moreover, mean 
square displacement and mobility correspond to the ex- 
ploration of the long time regime, while Cj and Ri decay 
quickly (a few collisions per particle) and therefore yield 
the short time behaviour. FDT could be valid in some 
time regimes and not in others. It must also be stressed 
that the measure of diffusion and mobility requires large 
times, so that averaging over many realizations becomes 
computationally very demanding and less precision is nu- 
merically available. In MD simulations in particular we 
have measured only the relation between C, and Ri. 



RESULTS 

In all the simulations performed (MD and DMSC) the 
values of the temperature ratios T1/T2 obtained are in 
good agreement with Q. This is not surprising for ho- 
mogeneous DSMC, while it is less obvious in MD (the 
agreement is mostly due to the low packing fraction n 
used) . 



concerned with steady-states, these extra terms are not 
present in the case studied. The other source of devi- 
ations, namely the deviation from Gibbs state, leads to 
small deviations 0] that could be difficult to detect . 



Particles 1 




• Particles 2 




y=x/(2T,) 




— -y=x/(2T 2 ) 










(a) 



oParticUs 1 




% Particles 2 




- v=.v/(27',J 




-y=x/(2T 2 ) 






(b) 



100 150 200 

B u (t) 



FIG. 1: Random force — 0), homogeneous DSMC: mobili- 
ties xi,2 vs. MSD -Bi,2; (a): an = 0.3, ai2 = 0.5, Q22 = 0.7, 
m 2 = 3mi, Ti « 0.2, T 2 « 0.38; (b): an = an = Q-22 = 0.9, 
m 2 = 5mi, Ti m 0.035, T2 fa 0.05. Symbols are numerical 
data, lines have slope l/(2Ti) and l/(2T a ). 
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FIG. 2: Random force (C = 0), DSMC with cells: mobilities 
Xi,2 vs. MSD Bi, 2 . (a): a w = 0.3, a 12 = 0.5, a 22 = 0.7, 
m 2 = 3mi, Ti « 0.1, T 2 « 0.185; (b): an = a 12 = a 22 = 0.9, 
m 2 = 5mi, Ti m 0.064, T 2 w 0.09. Symbols are numerical 
data, lines have slope l/(2Ti) and 1/(2T 2 ). 



Fluctuation-Dissipation ratio for the different 
components 

The measures of and Xi(t) allow us to check 

the linearity with t at large times: Bi(t) s» AD it and 
Xi(t) w fj.it. Moreover, figures H an d El clearly show that 
the Einstein relations are obeyed separately by the two 
components of the mixture, each with its own temper- 
ature: to the numerical accuracy one obtains x»W = 
■%pBi{t). Note that mobilities and diffusion coefficients 
of each species are different from their value in a monodis- 
perse case. We have considered various values of inelas- 
ticities, mass ratios, number density ratios, and kinds 
of thermostat, obtaining that this result is robust with 
respect to all these variations. We also note that a re- 
cent study [bgj has shown that the Einstein relation is 
not valid in its usual form in the case of an impurity im- 
mersed in an homogeneously cooling granular: the main 
reason is that extra terms arise due to the evolution of 
the granular temperature with time. Since we are here 



We now turn to the measure of R\(t) and Cj(t). Tech- 
nical details of the numerical procedure to perform this 
measure are given in [14|, From the definitions of £, Cj 
and Ri, it is clear that i?i(0) = l/(2m,*), and Ci(0) = 
= Ti/rrii. On the other hand, \im t ^oo Ri{t) — 
lirnt_>oo Ci(t) — 0. Thus we have plotted in figure EH 
the functions 2rrnRi(t) vs. niiCi(t). The FD relation 
Ri(t) = ^d(t) is verified replacing T by the partial 
granular temperature T\ and T2 of each component. The 
same result holds for different (rik in order to 
satisfy the properties of the perturbing force ©). 

The MD simulations (whose results are reported in 
figure 0) are more realistic since they include excluded 
volume effects and collision-induced correlations (which 
break the Molecular Chaos hypothesis, see for exam- 
ple However, at not too high packing fractions, 
we still observe the same results for FD relations. Larger 
packing fractions lead to strong heterogeneities in both 
density and granular temperature, giving rise to devia- 
tions from FD . 
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FIG. 3: Random force (i.e. £ = 0), DSMC with cells, and MD: 
2miRi(t)vs. m i C l (t) (i = l,2). Here n fc = 4. d{0) =Ti/rm, 
Ri(0) = l/(2m.j). In all cases the circles correspond to the 
FD plot for the type 1 particles, the pluses to the FD plot for 
type 2 particles. Straight lines have slope 1/Ti, dashed lines 
have slope I/T2. Top: DSMC simulations, (a): an = 0.3, 
a i2 = 0.5, a 22 = 0.7, m 2 = 3toi, Ti « 0.09, T 2 « 0.18; (b): 
an = ai2 = 022 = 0.7, 7712 = 5mi, Ti ft! 0.16, T2 « 0.35; 
Bottom: MD simulations, with n = 0.1. (c): Ni = N2, 
an = 0.7, ai2 = 0.8, a 22 = 0.9, m 2 = 3rm, Tj w 0.75, 
T 2 w 1.22; (d): iVi = 9iV 2 ; an = «i2 = Q22 = 0.9, m 2 = 5mi, 
Ti « 0.93, T 2 w 1.69. 

Fluctuation-Dissipation ratio for the whole system 

An interesting question concerns what happens to 
the Fluctuation-Dissipation ratio when measured for the 
whole system and not separately for the different com- 
ponents of the mixture. In other words one could ask 
whether one can define an effective temperature for the 
whole system and what is the relation of this temper- 
ature with the temperatures defined above for the two 
components of the mixture, or with the global tempera- 
ture T = X1T1 + x 2 T 2 (where x 1<2 = %^). 

A global measure would give, for the mean-square dis- 
placement 

B(t) = x 1 B 1 {t)+x 2 B 2 {t) , (7) 
and for the response 

X (t) = xixi(t) + x 2 X2(t) ■ (8) 



equipartition is satisfied. In figure^we present evidences 
supporting this view. 




B(t), BJt) 

FIG. 4: Homogeneous DSMC, random force, an = 0.3, ar2 = 
0.5, a 22 = 0.7, 7712 = 3mi, Ni = N-2. The filled squares 
correspond to the numerical data of x(t) vs. B(t). Circles 
and stars correspond respectively to particles of type 1 and 
2, and the thin lines have slopes l/(2Ti) and 1/(2T 2 ). The 
thick continuous line has slope 1 / (2T) while the dashed line 
has slope /V(4D) with 2D/n given by eq.H Here T ps 0.29 
while 2D /n w 0.27. 

Another way to consider the problem is to look at 
the ratio between R(t) = x\Ri(t) + x-2R2(t) and C(t) = 
x\Ci(t) + x^C^it) ■ In this case the impossibility to define 
an effective temperature reflects itself in the non-constant 
ratio between C(t) and R(t), as we show in figure |3J 




C(t), C.(t) 



Using the previously checked result that Xi(t) ~ A 4 ** 
and Bi(t) w AD{t : one obtains that B(t) w ADt and 
x{t) rs [Mt, with D = x\D\ + X2D2 and /i = Xi/j,i + X2[±2- 
The Einstein relation for the global case thus reads: 

2D = [X1D1+ x 2 D2)T 1 T 2 

(X X1D1T2 + X2D2T1 

It is clear that the ratio — corresponds to the global 
granular temperature only when T\ = T2, i.e. when 



FIG. 5: DSMC with cells, random force, an = 012 = a 2 2 = 
0.7, m 2 = 5mi, Ni = N2. The filled squares correspond to the 
numerical data of R(t) vs. C(t). Circles and stars correspond 
respectively to Ri(t) vs. Cr(t) and R2(t) vs. C 2 (f). Thin 
lines have slopes l/(2Tr) and 1/(2T2). The thick continuous 
line has slope 1/(2T); the dashed line is a guide to the eye 
showing that R(t) vs. C(t) is not a straight line. 

It should be remarked how these properties of the 
global Fluctuation-Dissipation ratio allow us to make the 
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FIG. 6: Circles: Mobility vs. Diffusion of a single particle of 
mass mtracer in contact with N = 500 particles of mass m, 
immersed in a heat bath (i.e. random kicks plus viscosity). 
We use the following conventions: atracer = Q12 and a = an. 
Only in case A) the tracer is also in contact with the external 
driving heat bath. A) m tra cer = m, a — 0.9, atracer = 

T tracer = Q yg. B ) mtracpr — 777,, OL = 0.9, 
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following prediction. Suppose to perform a global mea- 
surement of the Fluctuation-Dissipation ratio on a sys- 
tem of unknown composition. A global measurement on 
a monodisperse system would yield a well defined FD ra- 
tio and a well defined temperature equal to the granular 
temperature, independently of the observable used (see 
[l4j). On the other hand a polydisperse system would 
feature a typically non-constant FD ratio, function also 
of the observable (unless the system is completely elas- 
tic). 



The single tracer case 

Finally we investigate the special case N2 = 1, i.e. 
the case of a single tracer, immersed in a granular gas 
of N — Ni different particles, acting as a non-perturbing 
thermometer. FigureElreports the corresponding results. 
It turns out that, due to the inelasticity of collisions 
the tracer reaches a granular temperature different from 
that of the surrounding particles. Measuring a FD rela- 
tion thus yields the granular temperature of the tracer 
but not that of the surrounding granular gas. The tracer 
is thus sensing a temperature whose value is the out- 
come of the complex interaction between the tracer it- 
self and the granular gas. It is worth stressing how this 
happens even if the tracer almost does not perturb the 
granular. This leads to the non-conventional result that 
the "temperature" measured as a fluctuation-dissipation 
ratio |l5| . in a driven granular gas, depends on the ther- 
mometer. 



CONCLUSIONS 

In this paper, we have shown by numerical simula- 
tions that, in a binary granular gas, each component 
of the mixture obeys a FD relation with its own gran- 
ular temperature, while a global measurement provides 
a suitable definition of temperature only when equipar- 
tition is satisfied. These FD relations can be measured 
by different correlation-response pairs. In particular, the 
case of a tracer, which can act as a "thermometer" since 
it does not perturb the granular, has been investigated, 
and leads to the unusual conclusion that the measure of 
the temperature through FD relations would depend on 
the interaction between the tracer and the granular gas. 
Notice that theoretical studies (see e.g. (if) may then 
allow for an estimate of the gas temperature, knowing 
the temperature measured by the tracer. Further inves- 
tigations are needed to explore the scenarios bringing to 
violations of the FD relations [l^. This point becomes 
particularly important when, with more realistic energy 
injection through boundaries, we could expect hetero- 
geneities giving rise to position-dependent temperatures. 
It would be interesting to investigate how a tracer parti- 
cle would sample this inhomogeneous gas and ask about 
the meaning of the corresponding measured temperature. 
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